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ABSTRACT 



A general description of linear, swit ched-parameter 
systems is developed in terms of the mathematics of state 
variables. The restrictions imposed by the grouping of 
state variables and by' the "core" states of the system upon 
the coefficient matrices are considered. Two procedures 
for determining the element values of those matrices are 
each illustrated with an example. A general expression 
for a Cost function to be used to measure system quality 
is developed and illustrated with two examples. Extensive 
recommendations for future work are made. Several examples 
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for the improvement of the step responses of some switched- 
parameter systems are presented. 
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I. INTRODUCTION 



The first attempt to describe the behavior of a physical 
(real-world) system is usually based on a linear model of 
the system because such models are mathematically tractable 
and haVe been extensively studied. One purpose of this 
thesis is to develop a linear model using state variable 
theory for some real-world systems which may be modeled as 
linear systems with switched parameters. The parameters to 
be switched are the coefficients of the variables in the 
state variable representation of the system and the values 
the parameters take reflect the actions of various switches 
in the real-v;oria sysuem. Ailhuugh the model that is 
developed may be applied to some kinds of electronic cir- 
cuits, this thesis deals primarily xvith linear, swltched- 
parameter models of control systems. 

One example of the kind of real-world systems for which 
a model was sought is the arm positioning servomechanism in 
a magnetic disk-file. This servomechanism has two tasks. 

The first is to move the arm so that the read/write heads 
on the' arm move from one "track" to some other in the 
shortest time possible. A "track" is defined as being one 
of the narrow concentric rings on the disk's surface which 
is used to store information. The second task of the 
servomechanism is to maintain the read/write heads (hence, 
the arm) in a nearly constant position with respect to the 
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track of Interest. The arm Is positioned with respect to a 
reference set of pre-recorded tracks which reside on a 
special disk. 

During the transfer of the arm from one track to another, 
the servomechanism exhibits two modes of operation. While 
in the first mode, the system behaves as a "Bang-Bang" con- 
trol system, undergoing first a period of maximum acceler- 
ation, followed by a period of maximum deceleration. During 
the time that this mode is in effect, the position of the 
arm Is reckoned by counting tracks as they are passed. At 
the end of the deceleration phase of the first mode, the 
arm Is positioned In the near vicinity of the new track. 

At that time the servomechanism Is sivltched to the second 
mode wherein It benaves as a ilntai.-, high-gain, ccr;:pens at »o , 
closed-loop position control system. This second mode Is 
necessary in order to obtain the required tracking accuracy. 
The transition between modes Is accomplished by switching 
into the system a linear error detector and a compensator, 
and by disabling the "Bang-Bang" circuitry. 

The behavior of a system such as that just described 
may be modeled during each mode of operation by considering 
the "Bang-Bang" part and the linear • position control servo 
part separately. However, the behavior of the system during 
the transition from one mode to the next Is not descrlbable 
by either of the tv/o kinds of system alone. This thesis 
offers a solution to the problem by developing a general 
linear, sv/ltched-parameter model v/hlch Is applicable to 
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these kinds of systems and which allows one to deal much 

more directly with the transition between modes. The model 

developed uses state variable concepts and is reached as 

t h 

the limiting case of a general, n^ order, linear, time- 
varying system. The model is shown to be applicable to a 
restricted class of nonlinear systems. 

Previous work dealing with switched parameters in linear 
control systems has dealt primarily with the use of discon- 
tinuous damping factors in various kinds of servomechanisms 
[1-M. A different and particularly interesting application 
of swltched-parameter techniques is that reported by Flugge- 
Lots and Taylor [51 . The nonlinear system that they 
describe is a linear control system with multiple feedback 
paths and a relay switching mecnanism fui- determining v.’hich 
one of the feedback paths is to be used at a particular time. 
The work presented in this thesis is believed to be the first 
attempt to develop a model applicable to the above kinds of 
systems . 

The second purpose of this thesis was to develop a 
flexible and general means of measuring the "quality" of 
a particular swltched-parameter system. As an example of 
what is meant by the "quality" of such a system (and as 
another example of what is meant by "swltched-parameter 
system") consider the swltched-parameter version of a velo- 
city control servomechanism which is to be used to bring a 
shaft from a stopped condition up to a certain specified 
rotational velocity. The shaft is to be brought up to 
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speed by following a pre-planned acceleration curve until 
the shaft velocity is within the lock-in range of a linear 
phase-lock system which is to maintain the shaft at the 
desired velocity to within close tolerances. The inputs 
to the phase-lock system are a reference frequency and the 
output of a rotational-velocity-to-frequency transducer 
which is coupled to the shaft. At the Instant when the 
rotational velocity of the shaft reaches a certain value 
within the lock-in range, the system is to be switched from 
its acceleration mode to the velocity control servo mode. 

In such a system, one is concerned with how closely the 
programmed acceleration curve is followed, how much varia- 
tion in the velocity exists while the system is in the 
velocity control servo mode, and the extent and duration nf 
the switching transients which may occur when the mode is 
shifted. The "quality" of the system is some measure of 
how closely the system response approaches the design goals 
in each case. 

In order to have a single measure of the system's quality, 
the concept of a Cost function was taken from Optimal Con- 
trol theory and applied to these systems. Besides providing 
one number as an Indication of the quality of a switched- 
parameter system, the use of a Cost function permits' one to 
utilize function minimization techniques in the design of 
the system. The use of such techniques was desired in order 
to "optimize" the performance of such systems. For instance. 
Cost function minimization could be applied to the velocity 
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servo just described in order to minimize the transient due 
to the switching from the acceleration mode to the velocity 
servo mode. Such a technique might also be applied to 
determine the optimum charge on one of the capacitors in 
the compensator which is switched into the arm positioning 
servomechanism in the disk-file in order to minimize the 
switching transient in that system. 

The usefulness of the Cost function minimization technique 
is demonstrated in the Appendix where the step function 
responses of several simple swlt ched-parameter systems are 
"optimized" by determining an "optimum" switching time which 
minimizes an applicable Cost function. 

Preceding this demonstration is the development of the 
model in Section ii w.nicn is followed, in Section III, by 
the consideration of some constraints the physical world 
imposes upon the model and then by the presentation of two 
methods of determining the model from the real-world system. 
The tv;o methods are illustrated by application to two 
different examples of the type of swlt ched-parameter systems 
which have acceleration modes followed by position control 
modes. The applicable state equations have been solved on 
a digital computer to illustrate the behavior of the systems. 

The Cost function is developed in Section IV and illus- 
trated with examples based on Cost versus switching time 
for the two examples presented in Section III. Section V 
is a presentation of a summary and som.e conclusions which 
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are reached based upon the work done and the examples used 
to Illustrate it. Section VI presents several suggestions 
and recommendations of topics for further study. 
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II. MODEL CONSTRUCTION 



The starting point for the construction of a mathemati- 
cal model to describe the linear, swit ched-parameter system 
was the general, linear, time-varying (LTV) system described 



as follows: 












x(t) 


= 


A(t)x(t) + B(t)u 

f\j f\j f\j 


(t ) 




(1) 


y (t ) 

'V 


= 


C(t)x(t) + D(t)u(t) 

^ '\j % 'Xy 






with initial 


conditions, x = x 
<\,o 


(0), 


and for te [0 ,T] ; 


and 


where 














is 


the n-vector of 


the 


system states , 




x(t ) 


is 


the time derivative 


of X ( t ) , 

'\j 




u(t) 


is 


the m-vector of 


control Inputs, 






is 


the r-vector of 


output variables. 




A(t) 


is 


an n X n matrix 


of 


coefficients , 




B(t) 


is 


an n X m matrix 


of 


coefficients , 




ct 


is 


an n X r matrix 


of 


coefficients, and 




D(t) 

'V, 


is 


an r X m matrix 


of 


coefficients . 





The determination of the A, B, C, and D matrices is a 

'V; ^ 'V/ ^ 

straightforward procedure which is explained in many texts 
and which will not be considered here. In fact, the deter- 
mination of the coefficient matrices from the physical 
system will be less of a problem than may be indicated by 
the time-varying nature of those matrices in (1) due to the 
restrictions which are soon to be placed upon them. 
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It is well known [6] that the time variation of the state 
vector which is the general solution to the LTV system in 
(1) is given by 

t _i 

x(t) = Z(t)x^ + X(t) / X(s) -^B(s)u(s)ds (2) 

O; Or O»0 0 ^ 'V/ 

where X(t) is an n x n matrix satisfying the following 

vector differential equation and where X(t) ^ represents 

% 

the matrix inverse of X(t). X(t) is determined by solving 

% 

the system of differential equations 
X(t) = A(t)X(t) 

*\f 'Xj 'Xj 

over the interval te[0,t], with initial conditions X(0) = I, 
where I is the n x n identity matrix. 

Let the control Interval [0,T] be divided into N sub- 
intervals of arbitrary length with the boundaries of the 
subintervals denoted by the sequence of times 

{t^j ^2’ •••> 



where t., - t = T, and t. > t. , for 0 < 1 < N. 

No’ 1 — 1-1 — — 

Now let it be required of A(t), B(t), C(t), and D(t) 

'X; 0 ^ ' Xj 

that they all be matrices of constants for t ^ + 6 <_ t ^ t^_j_^ - 6, 

for 6 «min(t.,-, - t.). l.e. the coefficient matrices are 
1+1 1 

constants except for- a small interval about the subinterval 
boundaries. If t^ + 6 is denoted as tt and - 6 is 

denoted as then the set of differential equations 

describing the state of the system in the subinterval 
te[tt, niay be written as 

x(t) = A(t^,t-^^)x(t) + B(tJ,t-^^)u(t) 






( 3 ) 




% 

0 






with initial conditions x = x(t?). The matrix A(tt,t..,) 
represents the constant value of the A matrix in the partlc- 
ular subinterval, with a similar interpretation applicable 
to B(tt,t. In the next subinterval, the A matrix value 

would be represented by A(tt_^^ ,t^_^ 2 ) > with a similar repre- 
sentation for the B matrix. 

'X, 

Since the A matrix in any particular subinterval will 
not, in general, have the same element values as the A 

'\j 

matrix in any other subinterval, and since the systems con- 
sidered are LTV, it is necessary that the A matrix in the 

'\j 

transition region betx'/een two subintervals be LTV. If the 
A and B matrices in the transition region are denoted by 
A'(t) and B'(t) respectively, then the counterpart of (3) 

'V f\j 

in the transition I'egion may be wi-IILku ao 



x(t) = A'(t)x(t) + B’(t)u(t) (4) 



v/lth Initial conditions x = x(t.) and valid for t£[t.,tt], 

Oj 1 11 

The form of the solution to this set of differential equa- 
tions is given by (2). The value of the states at t = t? 
is of particular Interest since x(tt) serves as the initial 

Oi 1 

condition vector for the differential equations (3) which 
are applicable in the following subinterval, te [tt , tT_j_^] . 
Then from (2) the expression for x(tt) may be written as 



x(tt) 

'V 1 



X(t . )x(t . ) 

'X. 1 1 



xct;!") 

1 



/ 

t' 



X(s) 



-1 



B' (s )u(s )ds 

( 5 ) 



v;here, as before, X(t) satisfies 

O; 



X(t) = A' (t)X(t) 
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m 




4 



over the interval [t.,tt] with initial conditions X(t.) 

11 f\j ± 

and the end point values of A'(t) constrained such that 



I 



A' (tj^) - 

A'(t+) = A(t.,t^^p. 

Similar requirements may be imposed upon the B' matrix. 

Since the Interval [tT,tt] is a transition period during 
which the system coefficient matrices A and B alter value. 



'\j 



'\j 



it is interesting to determine how the value of the state 
vector x(tt) varies as the transition period is made smaller 

'Xt 1 

and smaller. In particular, what is lim x(tt)? Prom (5) 

6 -> 0 ^ 

it is seen that the Interval over which the solution of the 
X(t) equation is required is just [t7,tt]. This Interval 

^\j 11 

goes to zero length as 6 0. Hence 

11m X(tt) = X(tT) = I. 

6 -0 ^ ^ 'v 1 



Consequent ly , 






lim x(tt) = lim {X(tt )x(t . ) + X(tT) / X(s) ^B' (s )u(s )ds} 

u • 

1 



Or 



% 






Using the limiting values of tt and t. and XCtt), gives 

1 1 'X; 1 



ti 

lim x(tt) = Ix(tT) + 1 f X(s) ^B'(s)u(s)ds 

'Xj 1 'Xj'Xj 1 t . ^ ^ 



6-^0 



= x(t ) , 

^ 1 



( 6 ) 



since the value of the Integral taken over zero Interval is 
zero for a finite integrand. The Integrand will be finite 
so long as u(t.) is finite since B'(t) is assumed to be 
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finite and since X(t) ^ I as t7 ->■ tt. Hence, if the 
transitions between constant coefficient matrices are assumed 
to be no-elapsed-tlme switching of the coefficient values 
(which is accomplished by no-elapsed-time switching of 
system parameters) the values of the state variables will 
not be affected by the discontinuous change of the coeffi- 
cient matrices. 

So far just the particular problem of the transition 
between adjacent subintervals within the overall control 
interval has been considered. Hovjever, the limit expressed 
in (6) now permits the entire problem to be treated in a 
straight-forward fashion. Because the differential equations 
(3) for each subinterval are known, the solution of (1) 

(with the plecev;lse constant limitations on A, E, C, anci '.'j 

'X; '\j 

may be obtained by the sequential solution of N systems of 
n linear differential equations with constant coefficients. 
Thus, the differential equation 

x(t) = A(t ,t.^)x(t) + B(t ,t,)u(t) 

would be solved over the Interval te[t^,t^] with the initial 
conditions x = x(t ) = x(0). The value of the state vector 

'X^O a. O '\j 

at t = t, , x(t-), would then be used as the initial condl- 

1 'V 1 ’ 

tlon for the system of differential equations 
x(t) = A(t,,t,)x(t) + B(t ■,t,)u(t) 

which would be solved over the interval te[t^,t 2 ]- The 
final value of the state vector in this Interval, x(t^) 
would be used as the Initial condition vector for the next 
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subinterval. This process would then be repeated until the 
solution for the entire control Interval had been obtained. 

It has been assumed, to this point, that the coefficient 
matrices are all altered at the same Instant, at each of 
the subinterval boundaries. However, It Is clear from (6) 
that the limit will be unchanged even If the B matrix Is 
not altered or If It Is altered at some time either before 
or after the A matrix undergoes Its transition. The case 

aJ 

for which the A and B matrices undergo separate transitions 

rJ rJ 

Is conveniently viewed as a modification of the control 
Input, u(t), caused by the variation of the B matrix, applied 

/V 

to a system for which the A matrix elements are switched at 
the times t^^, •••> 

Tne system may of uuux'se be treated as other than 
sequential, linear, constant systems by using some partic- 
ular expression for the time variation of the A' matrix 

rJ 

during the transition period [tT,tt] subject to the con- 
straints on the endpoint values of the A' matrix mentioned 

A/ 

previously. It Is clear, however, that If the transition 
period Is a small fraction of the shortest natural time 
constant of the system, then the Instantaneously switched 
model of the transition will describe the behavior very, 
well. If the transition period Is not small with respect 
to the shortest natural period of the system, then the LTV 
problem must be solved for the transition period In order 
to obtain the Initial conditions for the next subinterval 
In which the system has constant coefficients. The speed 
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, 'i' 




of present day switching devices makes the no-elapsed-tlme 
transition a practical assumption for many systems. As a 
result, in the material that follows it is assumed that 
the transitions in all of the coefficient matrices are 
caused by an ideal (instantaneous) switching action. 

For completeness sake, it should be mentioned that the 
continuity of the state variables at the transition instants 
does not imply that the components of the output vector, y, 

rJ 

are continuous. Indeed, this is obvious from (1), since 
for x(t) and u(t) both continuous functions of time, any 
discontinuity in the C or D matrices will appear as a 
discontinuity in y(t). Let 






= + - 



C(t . T 

/vj T - 1 



,t, ) 



and 



AD. = D(t . ,t . . T ) 

rvl rJ 1 1+1 



D(ti_i,ti) 



rJ 



be the difference matrices for the switching at t^, where 
1 £ 1 <_ N-1. Then the vector Ay. = y(tt) - y(t7) will be 
the vector discontinuity caused by the discontinuous changes 
in the C and D matrices. Thus, 

aj /V 

Ay = C(t . ,t . , )x(t^) + D(t . ,t . , )u(t^) 

/\j 1 1+1 1 1 1+1 1 

- C(t. T,t.)x(t7) - D(t. .,t.)u(t7). 

rJ 1 — 1 1 1 pJ 1 — 1 1 1 

Having already assumed that x(t) and u(t) are continuous 
functions of time and, thus, do not change value at the 
switching instant. Ay. becomes 



Ay. = AC. x(t. ) + AD.u(t.). 



N- 



/x/1 ^ 



rJ ± N 



( 8 ) 
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Ifj however, u(t) Is allowed to change discontlnuously 
(x(t) is still constrained to be continuous) at the instant 
t^, with 



Au 



u(t^) 

N 1 



u(t ) 
1 



3 



+ + "* — 

then, substituting for u(t.) in (7), and with x(t.) = x(t.), 

± 1 /\^ 1 

gives 



AYi 



AC. x(t . ) + D(t ,t ., ) 

/V 1 ^ 1 1 It 1 



u(t . ) + Au. 

/N-' 1 1 



- D ( t . , , t . ) u ( t . ) 
1 — 1 1 ^ 1 



Hence, the difference vector becomes 




AC. x(t . ) 

1 /V 1 



AD. u(t . ) 

Af 1.AJ 1 



+ D(t^,t^^^) Au^ 



( 9 ) 



^ vr> v»/^ er> u 
^ 



( Q ^ 

s ^ * 



t* r.-p 



ni sr'.nm.i nuliy Oi .v v t ) 

rJ 



X j. v.^ 



general in the sense that it applies to any discontinuous 
change of the C or D matrices or the control vector u(t) 

rJ 



independent of any changes in either the A or B matrices. 

/V/ 

A change in these latter two matrices will have an indirect 
effect on y(t) since their change can be expected to 
influence the behavior of the state vector, x(t). 

/N/ 

The preceding work is founded upon knowing the answers 
to two questions. Both of the questions exist at the inter- 
face between the analysis and the design of the types of 
systems under consideration. That is, one must have first 
solved the synthesis or design problem before the two 
questions may be answered and before the analysis may 
proceed. 
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The first question arises because the preceding develop- 
ment assumed a priori knowledge of the switching times 
■tfi ^23 •••3 Because the manner in which these times 

are selected can completely alter the character of the 
system, the effects of the selection methods should be con- 
sidered. If the N -1 switching Instants are selected In 
some fashion prior to the beginning of the control Interval, 
then the analysis of the system conceived as a sequence 
of constant, linear systems would proceed as described 
above. The examples presented in Sections III and IV and 
Appendix A are all of this type. On the other hand, the 
switching instants may not all be fixed prior to the begin- 
ning of the control interval, but may instead be chosen 

^ ( r-s 4— ■» -r ■4— /— N ^ r>-4~ r-3 -? yv\ r-s /-3 3 1 -v^ -? -r- t-» r-s ^ 4r- s *1 -? ^ T T ”1 

with their selection being based upon some criteria involv- 
ing the current and projected state of the system. The 
real-world systems described in the Introduction are all 
examples of systems which determine the appropriate switch- 
ing Instant by comparing the system’s current state with a 
predetermined state which is fixed by the system designer. 
These systems are intrinsically nonlinear. The nonlinearity 
occurs since the functional dependence of the switching 
instants must be displayed as 

t . = t . (x(t ) , u(t ) ,t) . 

1 1 /V A> 

This functional dependence of the switching instants requires 
that the coefficient matrices be considered, in general. 
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not only functions of timej but functions of the state and 
control vectors as well. 

The functional dependence of the A and B matrices may 

/'J 

be displayed as : 

A = A(x(t) ,u(t) ,t ) 

B = B(x(t) ,u(t ) ,t ) 

/V/ 

with similar expressions for C and D. That the system Is 

/V /SJ 

Inherently nonlinear Is evident since the character of the 
system's response will be. In general, dependent upon the 
Input. However, the matrices are still constrained to be 
piecewise constant over the various subintervals of the 
control Interval [0,T]. Consequently, the general functional 
a-t-T An<3hi ps displayed above describe only formally the 
functional dependence of the very narrow class of systems 
which are considered here. 

If, once the switching Instant Is determined, the 
transitions of the coefficient matrices proceed without regard 
for the state or control vectors, then the limiting argument 
used In obtaining (6) shows that the states are still con- 
tinuous at the switching Instant, and that the linear, 
piecewise constant model can still be used to give precise 
results for what Is technically a nonlinear system. The 
applicability of the model to a variety of real systems Is 
thus assured. 

The second question that requires an answer before there 
can be a completed mathem.atlcal model probably has no 
"right" answer. It requires that one have determined the 
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best structure for the system which Is to be employed. If 
the control Interval is divided into two or more non-zero 
subintervals, then it is necessary to determine which 
sequence of coefficient matrices is the best. Alternatively 
one may seek that sequence of linear, piecewise-const ant- 
parameter control systems which would best accomplish the 
required task. This problem may be formalized as that of 
finding, for Instance, the A matrix of a particular system; 
l.e. find the function, f, such that 



A(t . ,t 

rJ 1 



1+1 



f (X(t) ,U(t ) ,t , A(t . n,t.)) 

rJ /\j /\> 
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for 1 < i< N-1, will cause the resultant system to have 
the desired characteristics. This problem remains to be 
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some sort of heuristic approach which is aided by the 
physical constraints that a particular application Imposes 
will lead to a system which gives acceptable performance. 

From the preceding it has been seen that linear, 
switched-parameter systems are easily handled by the mathe- 
matics of the state variable theory. It is apparent, how- 
ever, that the set of coupled differential equations that 
describes one of these systems is probably more amenable 
to being solved on a digital computer than to being expressed 
in a mathematically closed form. The closed form expression 
is possible in principle but would become Increasingly 
clumsy as the system order and the number of subintervals 
increased. 
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Another possibility is that hybrid computer techniques 
might be used to obtain a solution to the state equations. 
In applications that do not require a great deal of 
accuracy in the solution, such an approach may prove useful 
however, the insensitivity of the hybrid computer to small 
but significant changes in a parameter value or the length 
of a control subinterval may not allow the use of these 
techniques . 
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Ill, APPLICATION TO PARTICULAR LINEAR SYSTEMS 



A. THE RESTRICTED STRUCTURE OP THE SYSTEM 

The general swltched-parameter system which was the 
concern of the previous section will seldom be encountered 
In Its full generality. One reason for this Is that physi- 
cal systems exhibit a part of their reality as restrictions 
on the forms that the coefficient matrices of the state 
variable representation are free to take. As "an example of 
this sort of restriction may be cited that shown by the 
classical linear description of a servomotor. The transfer 
function Is shown with a block diagram of the servomotor in 
Flgui-e 1. 



u(t) 





1 






s ( s ♦ a ) 





y(t) 



Figure 1. Servomotor Block Diagram. 



The second order differential equation which describes the 
behavior of this servomotor Is 

y(t) + ay(t) = ku(t), (11) 

where the dots Indicate the time derivatives and a and k 
are constants. 

In putting (11) Into state variable form It Is common 
practice to choose as the states of the system the function 
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y(t) and its time derivative y(t). So let 



x-j^(t) = y(t) 

and 

X 2 (t) = y(t) 

be the system states. Then, rewriting (11) in state vari- 
able form: 



x^(t) = x^Ct) 

X 2 (t) = -ax 2 (t) + ku(t)j 



which may be put into matrix form as 



x^(t) 

X2(t) 



0 1 




x^(t) 


0 -a 




_X2 ( t ) 




( 12 ) 



(13) 



or 

x(t) = Ax(t) + Bu(t) . 

/\j rJ rJ 

The point of this elementary exercise has been to show that 
the forms of the A and the B matrices are restricted by the 
particular device and by the quantities selected as the 
state variables. In (13) the elements of the first column 
of the A matrix will always be zero for this particular 
choice of states to describe this particular device. Addi- 
tionally, the element of the B matrix that is presently 

zero, will always be zero under these conditions. Indeed, 

t h. 

one may say that for an n^ order system which has this 
servomotor in it and which uses y(t) and y(t) as the two 
states that describe the behavior of the servomotor (and 
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n-2 other states to describe the rest of the system) , If 
X. (t) corresponds to y(t) and x.(t) corresponds to y(t) 

(i 7 ^ j ) , then 

X. (t ) = X (t ) , (IM) 

l.e. the only non-zero term in the i^*^ row of the system's 
A and B matrices will be a. . = 1. The restriction on the 

rJ IJ 

elements of the two coefficient matrices result because the 
servomotor is described by a second order differential 
equation which, when expressed as coupled first order dif- 
ferential equations will always result in (1^) as the first 
order differential equation for x^(t) (= y(t)). Furthermore, 

the two states x.(t) and x.(t), or their equivalents, will 

^ <3 

always be present and always be paired as in (1^) in any 
system in which the servomotor is used (still assuming the 
basic y(t) and y(t) choice of states). 

B. THE CORE STATES 

The concept of grouping the system's state variables as 
occurred with x,(t) and x.(t) above forms the basis for the 
following discussion of swit ched-parameter control systems. 

No matter what the system's order may be, there will always 
■be a "core" of states which one is trying to control. For 
instance, in a shaft positioning servomechanism, one seeks to 
control the states of the servomotor which drives the shaft 
and those states (shaft position and shaft rotational velo- 
city) become the "core" states of the system. It may be 
that various compensation networks have been employed in 
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the system to improve its response. The state variables 
needed to describe the behavior of these networks, taken 
together with the "core" states, make up the n states 
necessary to describe the system. 

The essential point is that one seeks to control the 
values of only a few of the total number of the system's 
states. This does not mean that the remaining states are 
not to be controlled in some sense. It does mean that one 
is less concerned with their endpoint values than with the 
effect they may have on the variation of the "core" states. 

One is led naturally, then, to questions concerning the 
order of the variation of a particular state (or set of 
states) within a system. If the system can be represented 
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( 15 ) 



where, for an n^^ order system, A is an r x r matrix. A' 

A/O *J o 

is an r X (n-r) matrix. A, is an (n-r) x (n-r) matrix, and 
0 is a (n-r) x r null matrix, then, one may say that the 
states comprising the state vector x,(t) will exhibit varia- 

i. 

t h 

tions characteristic of an (n-r)^ order system since they 
are totally decoupled from the r states represented by the 
X (t) state vector. The states of the x (t) state vector 

A^O /\jO 

will exhibit the variations characteristic of an r^^ order 
system. 



X (t) = 



Ax ( t ) 

ajOaiO 



+ AiXj(t) + 






B u(t). 



28 



where the last two terms represent the forcing function for 
the Inhomogeneous system. 

The linear, swlt ched-parameter systems which are repre- 
sented In a particular time Interval by (15) are a subset 
of the more general class of swltched-parameter system 
described by 
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( 16 ) 



which Is the same as (15) except that the (n-r) x r null 
matrix Is generalized to the matrix A'. The x, (t) state 
vector for this case (as for (15) when that equation Is 
Interpreted as describing a swltched-parameter system) 
represents those of the system states which are considered 
the "core" states. 

The swltched-parameter system to be considered will be 
described by a sequence of N linear, constant-parameter 
systems during N successive subintervals of the overall 

t h 

control Interval which Is taken as [0,T], During the 1^ 
subinterval the system Is described by (l6) with the Initial 
conditions, x = x(tT) , and for te[t.,t.,T]. With A' and 

y O 1 1 1“T 1 1 

A' taken as non-null matrices and assuming that A Is rank n 

z^/O o A/ 

(non-degenerate system), the "core" state vector, x^(t) may 

have order behavior. It will be assumed that It does. 

The partitioning of the A matrix In (l6) serves to delimit 

the "core" states. The system's differential equations have 

been so ordered that such a partitioning Is possible. Now, 

assume that In the Interval [t..T,t..„] the behavior of the 

1 + 1 ^ 1+2 
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t h 

"core" states Is to be less than order. Such behavior 

could be caused by forcing the A' matrix to become [0;A, "] 

/\J± /V ‘ /vl 

s t 

In the (1+1) subinterval. This may require that the 
equations of the states not In the "core" be rearranged. 

The discussion will proceed as well If the A' matrix Is 
forced to become the null matrix. Thus, the behavior of 
the "core" states will be altered from n^^ order to (n-r)^^ 
order at t = It Is understood that the "forcing" 

of the A' matrix Is to be done by switching the values of 
the system’s parameters. The next linear system to follow 
that described by (l6), then, will be the system (15) with 
Initial conditions x = x(t7, , ) and valid for te [t , , , , t . . ^ . 

/v/ O /V X » X X » X X"T i:: 

The Initial conditions for (15) niay be viewed as the Initial 
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the Intervening sequence of linear systems to their present 
value, x(t. ,). The effect of the preceding systems Is, 

/vj X“r X 

thus, apparent. It may be noted that If knowledge of the 
variations of x (t) Is not required for further work, the 

rjO 

entire system of n differential equations need not be solved. 
One may solve for this special case just the equations 



x,(t) = A,x,(t) + B,u(t), 

/\/X /v»X/v/X a^Xa^ 



(17) 



with the Initial conditions being the appropriate subset of 
the Initial conditions for (15)- 

Since the order of the behavior of the "core" state 
vector Is (n-r), v/hlch Is assumed to be the number of the 
"core" states, the system v;lll be restricted In the changes 
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that may be made In it at the end of the current subinter- 
val. That is, the character of the behavior of the "core" 
states may remain of the same order, (n-r) , or the order may 
be increased (up to n^^ order) by suitably altering the A,' 
and the A' matrices to couple additional differential equa- 

rJ O 

tions to the equations of the "core" states. If the order 
is to be Increased, the values of those states in the x (t) 
vector which are to be coupled to the x^Ct) state vector 
at the next switching Instant will become quite Important 
since their values at the endpoint of the current subinter- 
val will be used as the initial conditions for the set of 
coupled equations which Include the "core" state equations 
in the next subinterval. Hence, one will be solving (17) 
augmented by some audlLional states, x'(t), ■■^here me 

a; O 

coupling is accomplished by the equivalent of the A^" matrix 

mentioned earlier as the non-null part of the A' matrix; 

1 

i.e. one must solve 

x’(t) =A,x, (t) + Al'x’(t) + B'u(t) (l8) 

-L ± ± ±r^ O ^ ^ 

where the initial conditions are given by ^.nd 

x'(t7,T) for the subinterval t e [t . . , ,t . . ^ , and where B’u(t) 

represents the appropriate forcing function terms from (15). 

It may be noted in passing that the initial conditions 

represented by x'(t7,,) may be controlled to some extent by 
Aj o 1 ' -L 

the system designer during the subinterval in which they are 
decoupled from the "core" states. This flexibility is 
allowed since the vector of states, x (t), in (15) inay be 

f\fO 
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treated as an Independent system so long as none of its 
member states are coupled to the "core" state vector, x^(t). 
The Improvements In the response of a swltched-parameter 
system that such flexibility may permit remain to be deter- 
mined. The values of x'(0) (for the case in which low order 

/V O 

behavior of the "core" states is desired during the first 
subinterval) are a particularly Important example of the 
flexibility that is available. 

Once the linear model of the system has been placed in 
the form of (15) > and the elements of the A matrix that are 
to change and the values that they are to take after each 
change have been established for each segment of the control 
Interval, then the solution of the system's differential 
ep u a u _L Or 1 .3 hi ay 

C. DETERMINING THE STATE EQUATIONS 

It is obvious from the foregoing material that the 
description of swltched-parameter systems is most easily 
accomplished using state variable concepts. Furthermore, 
if the state variables are restricted to those which have 
an immediate physical significance (e.g. capacitor voltages. 
Inductor currents, shaft rotational velocities, and the 
like) the conceptual understanding of the switchings of the 
system's parameters and the resultant behavior of the 
system's states is enhanced since one can readily see how 
much of one particular state is being coupled to some other. 
The order of the system and the values of many of the system 
parameters will have been fixed during the design of the 
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system. Usually the fixed parameters will be fixed because 
of the inherent physical characteristics of the system's 
components; although, some parameters may be fixed by 
"designer's choice." Those parameters that do not have 
fixed values during the entire control interval for one of 
the above reasons are those that are to be switched in value 
from subinterval to subinterval. 

Two procedures present themselves as possible methods 
for translating the physical system with swlt ched-parameters 
into the state equations of (16). The first procedure would 
have one determine the coefficient matrices (A, B, C, and D) 
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assuming a maximum amount of coupling between the state 
equations; that is, utilize individual multiplicative 
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the switching functions are non-zero, and determine a general 
set of coefficient matrices. Then, by allowing the appro- 
priate switching functions to take their zero values in the 
proper subintervals, one may read the element values of the 
coefficient matrices which are appropos of that subinterval. 

This procedure requires that one manipulate the differential 

t h 

equations for an n order system (assuming that the maximum 
order of the system is n), complicated by however many 
switching functions are required. For the case of N sub- 
intervals, there would be N-1 individual switching functions 
necessary . 

The second procedure asks that one proceed on a sub- 
interval-by-subinterval basis. Beginning in the first 
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sublntervalj the state equations of those states which are 
coupled to the "core" states are written down. Next, the 
equations for any other states which are not coupled to the 
"core" states but are coupled to a forcing function or 
other states not In the "core" (but are not Isolated) are 
written down, leaving only the Isolated states unaccounted 
for. One then considers the system as It Is to be after 
the first switching Instant, Identifies the Isolated states 
which have been coupled to the "core" states, and adds rows 
and columns to the appropriate system coefficient matrices 
In order to take these "new" states Into account. The 
elements of the coefficient matrices that correspond to the 
"old" states are changed appropriately so as to take account 
Kji uuupllii^ oi* L/iitr uncoupling u± Lliubo Cno cun — 

tlnues In this fashion, adding rows and columns (adding 
differential equations) to the coefficient matrices until 
they have attained their maximum dimensions, after which 
time the coefficient matrices are determined simply by 
altering the appropriate element values for each subinterval. 
Care must be taken to ascertain the Initial value of any 
"new" states when they are switched Into the system from 
their previous Isolation the first time. 

There are some pitfalls to be avoided In both of these 
procedures. For the most part, the pitfalls are related to 
the Identification of the states of the system. In switching 
parameters. It may occur that one obtains the equivalent of 
a circuit In which there Is a loop containing all capacitors. 
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or a node having only Inductors connected to it. Such a 
system is degenerate in that not all of the "states" repre- 
sented by the capacitor voltages, or the inductor currents, 
are Independent. This will require that the A matrix of 

rj 

the system in that subinterval have rank one less (assuming 
only one degenerate state) than it would have if the system 
were not degenerate. The system of equations may still be 
solved but methods that require inversion of the A matrix 

rsj 

must be used with care. The use of iterative differential 
equation solvers such as the Runge-Kutta and the Hamming's 
Modified Predictor-Corrector methods will provide time 
response solutions to the state equations with a minimum 
of concern about the values of the various coefficient 
matrices . 

Another type of pitfall may be encountered when using 
the first procedure mentioned above. The problem occurs 
when allov/ing a switching function to take a zero value 
causes an element of a coefficient matrix to become infinite. 
This is probably the result of dividing by a switching func- 
tion which was being treated as a symbolic, non-zero quan- 
tity during some algebraic manipulations. The problem is 
best cured by guarding against such manipulations. 

One other point needs to be mentioned. The system, in 
each sub-interval, appears as one or more sets of coupled, 
linear differential equations. In many cases the system 
being modeled villl be a classical, linear control system 
with a single input and a single output. Such systems are 
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well described using transfer functions, pole-zero plots, 
and the other well-known methods applicable to classical 
systems. Unfortunately, these classical analysis tech- 
niques are not well suited to the description of switched 
parameter systems even in the single-input, single-output 
case since the identity of the individual states is easily 
lost. This is not so important during the period that the 
system is in one particular subinterval, but it is very 
important at the ends of the subinterval because the values 
of the individual states must be known in order to determine 
the initial conditions for the subsequent subinterval. For 
instance, if in one subinterval a system has a transfer 
function, T(s) which changes as shown: 

T(s) = for t £ t^, 
s ^ 

and 

T(s) = — ^ 5 — for t > t., 

s (s+a) 

the behavior of the system's output will depend on the 
structure of the second system. If the system in the sub- 
interval prior to t^ is structured as three individual 1/s 
terms, then the system which results from the switching of 
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will have a different behavior than if the switching goes as 
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This will be so because the Initial condition on the l/(s+a) 
term will vary with its position in the system and thus cause 
the transient response term of the system output to vary. 
Since one must look at the individual state variables in 
order to determine the structure of the system in two adja- 
cent subintervals, it is convenient to describe the system 
only in terms of the state variables. Of course, it is 
still possible to consider the transfer function for a 
particular subinterval where that is possible and where it 
may lend additional insight into the behavior of the system 
in the subinterval. This will usually be advantageous when 
the effect of the parameter switching will be to bring in or 
to Isolate an entire group of states such as those descrlb- 

4^^ ^ ^ ^ ^ ^ ^ ^ 4- Vx 4 ^ 

u»xiv^ o Vw. J. '^-1_ kJ w I.XO o t.- CX X/ ky V., X XX 1 A Aq X k/XX_XkJ 

section . 

A simple example will be worked for each procedure to 
Illustrate its use. The two examples have been selected 
for mathematical simplicity and as illustrations of systems 
that show improved step response over purely linear systems. 
Both examples are third order systems with only one switch- 
ing of the parameters considered. Single-input, single- 
output systems were used with the output being one of the 
states in order to simplify matters even further and to 
remove the necessity for considering the C and D matrices 
in the examples. Because of the diagrammatic simplicity, 
the systems are shown in a block diagram form in the time 
domain v;ith each of the states represented by an integrator. 
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Figure 2. Block Diagram of a Third Order, Swit ched-Parameter 
Control System which has a Core of States Showing 
Third Order Behavior Before and Second Order 
Behavior After Switching. 

Figure 2 shows the system that is used to illustrate 
application of the first proceaure. Tne symool s^ represents 
the switching function and takes on the two values Indicated. 
The identity of the states, as shown, is arbitrary provided 
only that they remain distinct. The first procedure requires 
that all of the switching functions be considered non-zero. 
This is easily done by treating the functions (in this case, 
only s^) symbolically while writing down the n coupled 
differential equations that describe the system. Thus for 
the system of Figure 2, one writes 

x^(t) = (l-s^)u(t) 

x^Ct) = x^(t) 

x^Ct) = s^u(t)+ (l-s^)x^(t) + (l-s^)x^(t) 

- - s^x^Ct) 
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For equal to zero in the first subinterval and equal to 
one in the second, and with the switching taking place at 
t' = t^, the systems on either side of the switching point 
are readily seen to be 
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for (19) 
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x(t) = 


0 0 1 


x(t) + 
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u( t ) , for ( 20 ) 
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where the element values of the A and B matrices have been 



f\j 

developed explicitly. For this simple case the procedure 
works well. However, when the system complexity increases 
and there are more than several switching functions to con- 
sider, then the expressions required for the individual 
elements of the various matrices may become unduly unwleldly. 
The second procedure should then be used in preference to 
this one. 

The two sets of equations (19) and ( 20 ) were solved on 
a digital computer by means of an iterative routine that 
uses Hamming's Modified Predictor-Corrector method for 
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solving sets of linear differential equations. The method 
utilizes a fourth order Runge-Kutta procedure to generate 
four initial points, and then proceeds, using the preceding 
four points in an iterative procedure to generate an esti- 
mate of the next point, and, by extension, the solution to 
the set of equations. For further information on these 
methods see [7j 8]. Figure 3 shows the time response of 
the states of the system for two different values of the 
switching time, t^. The input, u(t) was a 10 volt step 
function in both cases. The length of the control interval 
was 10 seconds. The initial conditions of all the states 
were taken to be zero. The switching time for Figure 3(a) 
was t^ = 1.^5 seconds, and for Figure 3(b), t^ = 1.3^9 
seconas . 

Figure ^ shows the time response of the states of the 
second order system v;ith constant parameters, 

x(t) = 

to a 10 volt step input. Zero initial conditions and a 10 
second control interval apply also to this system. Figure 4 
is intended to serve as a reference for comparison with the 
time responses of the other systems. As may be seen from 
the comparison of Figures 3 and the switched-parameter 
system shov;s a marked Improvement in step response. 

The comparison of these two different systems is Justi- 
fied since the switched-parameter systems may be viewed as 
being the same second order system with means for modifying 
its initial conditions. Thus, one is comparing the response 
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Figure 3. Time Response of the Third Order, Swltched- 

Parameter System of the First Example to a 10 
Volt Step Input. (a) Switching Time, t^ = 1.^5 
Seconds. (b) t^ = 1.3^9 seconds. 
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Figure 4. Time Response of a Second Order System v/lth 
Constant Parameters to a 10 Volt Step Input. 
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of the constant-parameter, second order system when started 
from zero initial conditions to the response of the same 
system during the time its initial conditions are being 
"pre-charged" by the unstable third order system (in the 
first subinterval) and after it is released as a second 
order system (in the second subinterval). 

The choice of switching time, t^, for Figure 3(b) was 
decided by repeated trlal-and-error solutions of the state 
equations. The value obtained, t^ = 1.3^9 seconds, gave 
the best-looking response. An Increase or decrease in this 
value by as much as 0.001 seconds produced a visually 
detectable Increase in the overshoot or undershoot respec- 
tively. A similar sensitivity to the value of the switching 
tirut; riHH ucen noticed in all of L-hc examoles v'.".'±ke>l. TVie 
effects of such sensitivity may preclude the use of hybrid 
computer techniques in the solution of the state equations 
of such systems. 

The second example will be equally simple. Figure 5 
shows the third order system in the first subinterval. As 
before, the states have been individually identified and 
represented by Integrator blocks. 
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Figure 5 . Block Diagram of the Third Order, Swltched- 

Parameter System of the Second Example Prior to 
Sv;itching . 
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The system differential equations in this subinterval are 
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( 21 ) 



It is desired that the system take the configuration shown 
in Figure 6 in the subinterval te[t^,T]. The changes to 
be made can be systematically determined by examining the 
differential equation of each state in turn, as they are 
given in C21). 
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Figure 6. Block Diagram of the Third Order Switched- 

Parameter System of the Second Example After 
Switching. 



Going from top to bottom of the state vector, examining 
x^(t), X 2 (t), and x^(t) in turn, it is apparent that the 
new set of state equations becom.es 
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If an isolated fourth state had been encountered in the 
second subinterval , and if the coupling between the states 
x^(t), X 2 (t) and shown in (22), then the fourth 

order system of this subinterval would be formed as 
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where the elements of A,’ , At, A', and B, would be determined 
from the coupling between X|^(t) and the three prior states. 

It should be well noted, that since this is the first time 
that X|j(t) has been used in the system (else it would already 
be included in x(t)), that particular attention must be paid 
to the value of Xj^(t^), which is the initial condition of 
x^(t), for the new subinterval. 

Depending upon the method used to solve the ' set of 
differential equations (22) or (23), the chore of program- 
ming a digital computer for the task of solving those 
equations may be eased by always working with an n^*^ order 
system, even v;hen there are isolated states that are being 
"ignored." It would then be necessary only to alter the 






element values in the coefficient matrices and not necessary 
to change their dimensions during the course of the problem 
solution. The isolated states are conveniently allowed for 
even if not individually identified, by adding the equiva- 
lent of A '.A,, AJ , and B, as null matrices. If there is 
some trouble encountered in determining the maximum system 
order, the lov;er order system's equations may be written 
down and then adjoined to the appropriate null matrices, 
when the maximum system order is known. 

As for the first example system, this system, as des- 
cribed by (21) and (22) was solved on the digital computer. 
The trajectories of the states are shown in Figures 7(a) 
and (b) for a control interval of 10 seconds and a 10 volt 
step I’unccion input. Fl^uit; 7(b) is for a switching 
t^ = 1.3^9 seconds, and 7(a) is for t^ = 1.699 seconds. 
Figure 7(a), when compared with Figure 3(a), shows that 
the response of this system has a slightly longer rise time 
than the first system for about the same amount of over- 
shoot. The response of this system is also markedly better 
than that of the plain second order system. 
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Figure 7- Time Response of the Third Order, Swltched- 

Parameter System of the Second Example to a 10 
Volt Step Input. (a) Switching Time, t^ = 1.699 
Seconds. (b) t^ = 1.3^9 seconds. 
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IV. MEASURES OF SYSTEM QUALITY 



The various methods available for taking the measure of 
a system (or Its model) have the common goal of providing 
the analyst or designer with a quantitative estimate of hov; 
well the system may be expected to accomplish Its designed 
tasks. Linear system models are commonly tested by applying 
a known Input (usually a step function or a sinusoid) and 
comparing the output with the desired output. In the case 
of the step function Input, one might examine the rise time, 
overshoot, and settling time of the output response. As 
may be seen from Figures 3 and 7, there seems to be justi- 
fication for using rhose ineasui-es of sybLem quality to 
grade some swltched-parameter systems. 

Unfortunately, the usefulness of the classical system 
response measures Is somewhat degraded since It Is usually 
not possible to make the direct correlation between the 
classical response measure and a system parameter that Is 
possible with the classical linear system. The restriction 
results from what may be considered an Inherent ambiguity 
in linear, swltched-parameter systems. The ambiguity arises 
because there may be more than one sequence of linear, 
constant-parameter systems that would have the same output 
for a given input. This Is another manifestation of the 
intrinsic nonlinearity of sv/it ched-parameter systems. Still, 
for those swltched-parameter systems that are designed to 
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replace a classical linear system, the classical response 
measures may well provide a useful means of grading the 
systems . 

There are cases, however, for which the classical 
measures of system performance are not adequate. As one 
example may be cited the Linear Regulator problem of Optimal 
Control theory. A Cost function approach to the m.easure- 

ment of the quality of linear, swlt ched-parameter systems 

/ 

has several features to recommend it. The possibility of 
having several subintervals in which such systems may have 
greatly varied behavior and in each of which the designer 
may..-h.ave separate goals , strongly recommends a flexible 
measure of system quality. In addition, a Cost function 

bt; ubeTul In px'Oviuin^ a lucdbuic OT britVci'cLl VaricLulco 
in the system at the same time. For instance, one might 
construct a Cost function for a system whose value would 
be the combined measure of the system's settling time and 
the integral of the error between the output and the 
input squared (assuming a system for which such things 
make sense.) Another reason for desiring to use Cost 
functions with swltched-parameter systems is that such an 
approach lends Itself to some sort of system optimization 
accomplished through minimization of the Cost function. 
Several examples of Cost function minimization are presented 
in Appendix A as a demonstration of the utility of the 
method. The examples are limited to the determination of 
the optimum switching time for some simple systems but they 
illustrate a method that might profitably be employed. 
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Against this potential flexibility must be balanced the 
loss of detail that is inherent in the use of any sort of 
Cost function. The loss of detail arises, of course, 
since the Cost function is the weighted sum of several 
different quantities and there is always the possibility 
that an Increase in one quantity was offset by a decrease 
in another. Ambiguity in the value of the Cost function 
is possible when the changes in the different quantities 
offset each other so that the function value remains 
unchanged . 

A general Cost function for swlt ched-parameter systems 
may be written as 

+ RT^ (211) 

where 

Q = [Q]^,Q 2 ) • • and ^ 0 for 1 = 1,N, 

R > 0, 



are weighting factors for the Cost vector, J , and the 
length of the control interval, T , respectively. The 
vector 

is ^ • • ‘ 



the length of the control Interval is defined as 
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c 



N 

E 

1=1 



(t^-t 



1-1 



) 






for the sv/lt ched-parameter system given by 
x(t) = A(t)x(t) + B(t)u(t), 

/V/ />J /V 

and valid for t in the Interval [0,T^]. 
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In the 1^^ sub interval, 

t. 



A/ 



m 1+1 'V 'T’ 

J, = w(t . , ) H. w(t . , ) + / tv(s) M.v(s) + u(s) N.u(s)]ds 

1 1+1 /\J±^ 1+1 +. /V /'^l/^ 

^ (25) 



v(t) = x(t) - r. (t) , 

/v' 1 



where H. , M. , and N. are n x n, n x n, and m x m positive 

semi-definite diagonal matrices of weighting factors, 

respectively. The vector, r.(t), is the desired trajectory 

for x(t) in the subinterval te[t.,t..,], and x„(t.,,) is 
~ 1+1 ^ /^ D \ 1+1 

the desired value of x(t) at t ~ t.,,. An equo.tion analgous 
to (25) is applicable to each of the N subintervals. 

Which of the weighting factors are to be non-zero, and 
wnar une values are co oe is left to the system analyst 
for decision since the weighting factors should be chosen 
to provide a measure of those things that are of interest in 
a particular system vihen it is used for a particular purpose . 

Although a general statement of the useful forms of the 
Cost functions for switched-parameter systems has not yet 
been made, there are some indications that the full gener- 
ality of ( 2 ^) and (25) should be restricted depending upon 
the various applications. If the control interval is of 
fixed length for a particular problem, the weighting factor, 
R, may as well be zero since the RT term in (2^) will then 
only add a constant to J^. In a particular subinterval 

(the 1^^, say) the matrix M. would have non-zero elements 

1 

for each state whose deviation from a specified trajectory 
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was to be penalized (i.e. Increase the Cost function value). 

Llkev;lsej N. would have non-zero elements for those elements 
1 

of the control vector, u(t), whose scaled magnitude was to 

/V 

contribute to the cost In the subinterval. A similar state- 
ment applies to and the deviation of the values of the 
state variables at the endpoint of the subinterval from 
some specified goal, Xj^(t^_i_^). Of course. If the value or 
variation of a system variable Is to be Ignored, the cor- 
responding element In the appropriate weighting factor 
matrix Is set to zero. 

The Cost function measure of system quality may be used 
In .a function minimization scheme to optimize the variable 
parameters In the system. In so doing, one treats as a 
runctiori of the pc.r’O.inct ore to bo vsrlod ?r_d rridiilr^l^os 
function to obtain "optimum" values of the parameters. For 
example, one might optimize the sv;ltchlng Instants (the 
t^, for 1 = 1,N) for a particular swlt ched-parameter system 
with a particular Input, or one might try to minimize the 
total control Interval necessary for the system to complete 
Its response to another Input (by suitably varying the t^^, 
or the Initial condition due to an Isolated state, perhaps), 
or attempt to optimize one of the system parameters. 

In general, the possible values of .the Independent .vari- 
ables of the minimization will be limited (e.g. the t^ may 
be constrained such that 0 £ T, for 1 = 1,..,N) and 

the minimization of must use those procedures applicable 
to constrained systems. It may be noted. In passing, that 
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the control vector for the swltched-parameter system may 
be chosen so as to minimize a functional of the form of J,p. 
This represents the optimal control approach to determining 
the control function for swltched-parameter systems. It 
is an area that remains to be explored. 

As examples of their use, two Cost functions were con- 
structed for the system of (19) and (20) and Figure 2. The 
first Cost function was 



Jt 



Xi(Tc) - 0 

x^(T^) - 10 

X3(Tc) - 0 



0 

1 

0 



Xi(Tc - 0 

x^(T^) - 10 
X 3 ( T ^)-0 



( 26 ) 



Jip was treated as a function of the length of the first 

X. • 1 4- ^ ^ c? ^ ^ - A- , , 4 . _ 

OLt.k.>ULXAw«.....t. L/-I • vC.' VjJ- Ullt- t ( i L/ 

JU 

(19) and (20), over the control interval T^ = 10, and eval- 
uation of (26) for 100 incremental values of t^ spanning 
the Interval [0j2j , with a 10 volt step function Input to 
the system produced the points for the curve plotted in 
Figure 8. 

Figures 9 and 10 show the Cost function of (26) modified 
so that deviation from the desired final value during the 
entire second subinterval is counted as part of the cost. 

The modified Cost function is, then. 
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Figure 8. Cost Function of a Third Order, Sv/itched-Parameter 
System as a Function of the Time to the First 
Sv;itching for a 10 Volt Step Input and Where Only 
the Square of the Deviation of the States from 
the Desired Final Value is Counted as Cost. 
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Figure 9. The Cost Function of Figure 8 Modified to Count 
the Time Integral of the Square of the Deviation 
from the Desired Final State During the Entire 
Second Control Subinterval as Cost, for System 
Input Steps of Various Amplitudes. 
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Figure 10. Same as Figure 9 Except that Time to the First 
Switching is Allowed to Vary in the Interval 
[0,51 . 
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( 27 ) 



ds 



The scalar is both the desired trajectory (a constant) 
and the desired endpoint value for the second subinterval. 
It is also the amplitude of the input step function. 

The purpose of (26) and the first term of (27) is the 
same; namely, to penalize (Increase the value of the Cost 
ruiicLluii) olitr o±* tlid* sto-tcro Trcm z'lZZmT 'tcLr'gct 

values at the endpoint of the control interval. The Cost 
function, J^, is increased by the square of such deviation 
as does exist at that time. In both (26) and (27) the 
values of x^(T^) and x^(T^) are to be penalized if they 
differ from x^ or zero respectively. The endpoint value 
of x^(t) does not affect the value of the Cost function 
since the first element on the diagonal of the weighting 
factor matrix, , is zero. 

The second term of (27) serves to penalize the system 
if the trajectories of its states vary from some prescribed 
trajectory during (in this case) the second subinterval. 
Since the only nonzero element of corresponds to the 
deviation of x„(t) from x , , only the integral of the square 
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of that deviation will add to the Cost function. The result 
is that both overshoot and slow rise time of X 2 (t) in the 
second subinterval for the swltched-parameter system being 
considered will be penalized. 

If (27) were to be applied to a shaft position control 
system with constant parameters, and if X 2 (t) and x^(t) were 
the state variables corresponding to the shaft position angle 
and its time derivative respectively, then the application 
of the Cost function may be seen to give a measure of how 
closely the system approaches the value of the input step 
function after some period of time which is considered as 
the control interval. The Cost function would also give 
a measure of how rapidly, without producing excessive over- 
snoot 01 X/he xxri ct j- V cx i/iie V ' — fjL n a x. 

value. It may be seen from consideration of this example 
that the Cost function is ambiguous in that an overdamped 
system which rises slowly to its final value may have the 
sam.e value of cost as does an underdamped system which over- 
shoots and then oscillates about the final value several 
times . 

The multiple curves in Figures 9 and 10 show how the 
Cost function varies as the amplitude of the input step in 
varied. Figure 9 is analogous to Figure 8 in that t^ is 
only allowed to vary in the interval [ 0 , 2 ]. Figure 10 
shows hov7 rapidly the Cost function increases as t^ is 
increased beyond 2 and provides a non-rigorous basis for 
believing that the minimums evidenced by the functions in 
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Figure 9 truly global minimums . The system's states 

were all started from zero Initial conditions In the compu- 
tation of these curves. The possibility of minimizing the 
Cost functions of Figures 8 and 9 with respect to the 
"variable" t^ is apparent. 

The determination of Cost function plots such as Figures 
8, 9j and 10 requires a large expenditure of computer time 
since repeated solutions of the state equations are neces- 
sary to evaluate the Cost function for different values of 
its Independent variables. Each individual curve in the 
three figures required approximately three minutes of compu- 
tation time on the IBM 360/6? at the Naval Postgraduate 
School's W.R. Church Computer Center. These plots are for 
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would have been required for more than one independent 
variable. Fortunately, function minimization usually 
requires fewer evaluations of the Cost function and, hence, 
less computation time. 
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V. SUMMARY AND CONCLUSIONS 



The preceding material has 1 ) developed a general 
description of linear, swltched-parameter systems In terms 
of the mathematics of state variables, 2) presented the 
concepts of state variable grouping and of the "core" 
states of a swltched-parameter system, 3) presented two 
procedures for determining the sequence of coefficient 
matrices In the state variable representation that describes 
the behavior of the switched parameters In the system and 
Illustrated their use with simple examples, and 4) presented 
one possible Cost function which may be used to measure 
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examples to demonstrate Its application as a measure of 
system quality. That such a Cost function may provide a 
useful way to optimize the behavior of the system has been 
demonstrated In Appendix A. 

Several conclusions may be drawn from the work just 
summarized and the examples which were presented In conjunc- 
tion with that work. Although the model that was developed 
was not applied to a real-world system. It is obvious from 
both the mathematics and the examples that physical systems 
like those presented In the Introduction may. Indeed, be 
modeled as swltched-parameter systems. A conclusion that 
is possible as a by-product of the exam.ples of swltched- 
parameter systems that were used is that It Is possible to 
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greatly Improve the step response of a heavily damped second 
order system by using the states in the second order system 
as "core" states in an unstable second or third order 
system in order to "pre-charge" them. The penalty that one 
must pay for such improved response is a marked dependence 
of the response on the time that is taken to "pre-charge" 
the states. 

One may also conclude that it is possible to use Cost 
function minimization to, determine the optimum switching 
time which produces the "best" step function response in 
at least one variety of switched-parameter system. The 
optimum switching time so determined is dependent upon the 
form of the Cost function which was minimized. Finally, 



it lb evident that the- digital computer .methods used to 
solve the state equations and to determine the optimum 
switching point were both sensitive to and able to cope with 
the sensitivity of the system response to the switching 
time . 



This first attempt to study linear, switched-parameter 
systems in some generality has raised many more questions 
than it has ansv/ered; indeed, much the greater part of a 
complete study of these systems remains to be done. Not. 
all of the possible paths that lead from this tentative 
beginning have been discovered. Some of those which have 
been and which may prove v/orthy of further study are pre- 
sented in the following section. 
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VI. RECOMMENDATIONS FOR FURTHER WORK 



The areas of further study that are recommended below 
constitute a study of what may be called the "design prob- 
lem." Some of the questions require that one attempt to 
answer the question, "What is the optimum system to use for 
this job?" That question has been answered heuristically 
in the past and a definitive answer is not expected soon. 
Still, one may ask questions about the system architecture 
that do have some hope of being answered. 

The first problem relating to the architecture of the 
systems that might be studied would be to define the con- 
straints uii the physical system that are required to ohr.Min 
a certain order of behavior of the "core" states. For 
instance, one might ask what system architecture is required 
to produce third order behavior in a system v;ith two "core" 
states . 

The kind of behavior that one might seek to establish 
in a swit ched-parameter system in order to satisfy the design 
goals should be related to what is possible, desirable and 
obtainable in the way of system architecture. This might 
be approached by seeking to determine what sort of output 
responses are possible when one switches betv/een two par- 
ticular linear systems. If one were to form a table of 
systems before and after switching, and the particular 
response characteristics obtained, one might enter the system 



62 




I 




used in Figure 3 , noting that it goes from a third order 
system with a pole in the right half of the s-plane to a 
damped second order system and that the result is an im- 
provement in the system step response for a suitable 
selection of switching time. The hope is that [by looking 
at such a table] , one would be able to decide which system 
to switch to in order to obtain a desired response. 

A ’’building block" approach may also be possible. Prom 
the table just mentioned it may be possible to identify 
certain sequences of linear systems which have character- 
istics such that one might desire to use these subsequences 
as building blocks in constructing various systems. Of 
course, this raises the question of determining what sorts 
oj. building blocks iiu.g,nt be i'cmuxi'ccl , oi’ desareo, icr such 
a system and how to put them together to obtain the desired 
results . 

The work of Plugge-Lotz and Taylor [6] that was mentioned 
in the Introduction indicates that a quasl-adaptlve scheme 
can yield useful results. A generalization of those results 
in terms of swit ched-parameter systems would be desirable. 

For swit ched-parameter systems generally, one might 
inquire what the switching criteria should be. As an example 
one may not want to switch to a particular system while a 
capacitor has a non-zero charge. One might switch Instead 
(with the switching based on the values of some other states) 
to an Intermediate system which would discharge the capaci- 
tor before making the transition to the particular system 
first mentioned. 
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As a grand goal of the work with system architecture, 
one might seek a procedure for the systematic design of a 
switched-parameter system to meet some specified require- 
ments . 

Optimization problems represent another area of study 
that could produce significant results. One usually encoun- 
ters such problems in connection with optimizing an input 
to an already fixed system; however, with a system which is 
not fixed (e.g. switched-parameter) it may be reasonable 
to ask what sequence of systems gives an optimum response 
for a fixed input. This has the nature of an inside-out 
optimal control problem where one seeks to determine the 
system that makes a given control optimal. 

The pobsiuilil^ Of varying the initial conditions of 
certain of the states in the system was shown in Section III 
to allow a large Improvement in the response of a heavily 
damped second order system, and has been mentioned as a 
possible candidate for an optimization technique. This 
idea should be explored further, as should the theoretical 
Improvement in system response that may be attained. Being 
able to vary the initial condition of one of a system's 
states would seem to permit one to use the transient response 
in the system caused by that initial condition to alter the 
transient response of the other states of the system in a 
controlled fashion. 

The kind of Cost function that one uses will influence 
the results of the optimization of a system. The 
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correspondence between changes in the Cost function and 
changes in the physical response of the system is very 
important. As was seen in the various examples, a smooth, 
slowly varying Cost function may correspond to a physical 
response which is extremely sensitive to the particular 
parameter variation, as was the case with t^. The estab- 
lishment of guidelines for the assignment of weighting 
factor values, thus, seems useful. It may be worthwhile 
to investigate the type of Cost function necessary to allow 
a system parameter to be optimized with the result that 
some of the system's states may be constrained to stay 
within certain limits. This would be useful since it would 
allow the designer to specify physical limits on the system's 
states where such limits resuired from saturation liinlts on 
the physical device or from some desired performance criteria. 

It should be possible to formulate an optimal control 
problem for switched parameter systems. This should be 
determined and any limitations on the system imposed by the 
requirements for the existance and uniqueness of the optimal 
control systems with constant parameters and with the "opti- 
mal" control system mentioned previously in which the 
control 'is obtained by switching parameters. 

In describing the switching curves of the present relay 
control systems, it is common practice to speak of the 
switching boundary as a hyper-surface in state space. Fur- 
ther work might be done to determine the effects (beneficial 
or deleterious) of being able to sv^ltch on a particular 
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"line" in that hyper-surface, where such flexibility would 
be allowed by the Independent control of the initial condi- 
tion of a particular state which was isolated from the 
"core" states at some point. Such considerations could 
also be the basis for the design of certain small segments 
of a larger system. 



66 








k. 




APPENDIX A 



APPLICATION OF COST FUNCTION MINIMIZATION 



The purpose of this appendix is to demonstrate the 
utility of a Cost function minimization procedure as a 
means of optimizing the step function responses of some 
simple swlt ched-parameter systems. As the first example 
of a Cost function to be minimized, consider Figure 12. 

The set of functions shown there resulted from the appli- 
cation of five different amplitude step functions to the 
swit ched-parameter system of Figure 11, which is the system 
of Figure 2 reproduced here for convenience. The system 
was cjtditeJ from zero initial ccr.di'^icns for ?>T' in 

each case. 



u(t) 



^1 = 



te [0 




■X2(t) 



Figure 11. Replica of Figure 2. Block Diagram, of a Third 
Order, Switched-Parameter System. 
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Figure 12. Replica of Figure 9- The Cost Function of a 
Third Order, Swltched-Parameter System as a 
Function of the Time to the First Switching for 
Various Amplitude Step Inputs; where the Square 
of Both the Deviation of the States from the 
Desired Endpoint Values and the Deviation of 
the States from the Desired Final State During 
the Second Subinterval are Counted as Cost. 



68 



n n n 









Prom Figure 12 it is easily seen that the Cost functions 
each have a minimum for a time to first switching, t^, of 
about 1.35 seconds. For low amplitude steps (e.g. 5 volt, 
or 7 volt) the minimum is quite broad. It becomes more 
sharp as the step amplitude Increases to 20 volts. Because 
of the broad minimum, one would expect the time response 
of the system not to vary drastically for a small change in 
t^. A comparison of Figures (13(a) and (b) (which are 
replicas of Figures 3(a) and (b)) shows that this is the 
case. The small increase in the Cost function which is 
barely discernable on the 10 volt input step curve of 
Figure 12 as t^ increases to 1.45 seconds, is seen to be 
reflected in an overshoot-and-decay response in Figure 13(a). 

Not to bo tgnor*0'C!. ts tino I*2.Ct 
functions of Figure 12 all have about the same minimum. 

Hence, for this system, if the switching time is taken as 
1.35 seconds, the system response should be near optimum 
for a wide variation of step Inputs. The system behaves 
as a linear system, the character of its output response 
being unaffected by the amplitude of the input. This linear 
response is to be expected since the system is linear in 
both subintervals and the switching Instant, t^, is inde- 
pendent of the values of the states, being fixed in value. 
This is brought out vividly in Figures l4, 15, and 16 , which 
show the output response of this system for input steps of 
5, 7, and 15 volts respectively. 
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Figure 13- Replica of' Figure 3- Time Response of a Third 
Order, Sv;it ched-Parameter System to a 10 Volt 
Step Input. (a) Switching Time, t^ = 1.45 
Seconds. (b) t^ = 1.3^9 seconds. 
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Figure 1^. 



Time Response of the Third Order 
Parameter System of Figure 11 to 
Step with the Switching Time, t^ 
Seconds . 



Switched- 
a 5 Volt Input 
= 1.3^75 
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Figure 15. Time Response of the Third Order, Switched- 

Parameter System of Figure 11 to a 7 Volt Input 
Step v/ith the Switching Time, t^ = 1.3^75 
Seconds. 
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Figure l6. Time Response of the Third Order, Switched- 
Parameter System of Figure 11 to a 15 Volt 
Input Step with the Switching Time, t^ = 

1. 3^75 Seconds . 
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These figures should be compared with Figures 13(a) and 



(b) . Note that there is a slight change in the value of 
t^ from 1.3^9 in Figure 13(b) to 1.3^75 Tor Figures 1 ^ , 15, 
and l6 . 

The second example of the application of Cost function 
minimization will be for a second order system with a 
switched damping factor. The block diagram of the system 
before and after switching is shown in Figure 17. The 
state equations for this system are 



[u(t)] for 

te [0 ,t^] , 



x(t ) = 


'*0 


l” 


x(t) + 


" o " 




_0 


0 




1 



u(t) 




-x,(t) 



(a) 



te [0,t^] 




(b) te[t^,10] 



Figure 17. Block Diagram of a Second Order, Swltched- 
Parameter System. (a) Before Switching, 
(b) After Switching. 
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and 



x(t) = 


0 


1 


x(t) + 


"o" 


r>j 


_-l 


-i 


A/ 


1_ 



[u(t)] for 
te [t ^ 



10 ] , 



where the length of the control Interval Is taken as 10 
seconds . 

A Cost function for this system was constructed analgous 
to the one for the third order system just discussed, the 
function being dependent on the final values of the states 
at the end of the control Interval and dependent upon the 
deviation of the state x^(t) from the desired final value 
during the second subinterval. Thus, the Cost function was 



x^(10)-x^^ pi O^px^dO-x^j 

I j n , 

X2(10)-0j \_0 lJ[_X2(10)-0j 



10 



/ ( 




T 


r* — 1 


\ 




x^(s)-x^ 




o 

1 — I 


X^(s)-X^ 1 


l\ 


X2(s ) - 0 




_0 0 


x^(s) - 0 j 



ds 



( 28 ) 



The first term of (28) penalizes the deviation of the 
states from their desired values at the end of the control 
interval by increasing the Cost function by the sum of the 
squares of the deviations of both states from the desired 
end point values. Deviation of x^(t) from the input step 
amplitude (the desired final state for x^) represented by 
x^, and x^(t) from zero final amplitude are both penalized 
in this v/ay . The second term penalizes any deviation of 
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x^(t) from the desired trajectory (In this case, a constant 
value, x^, which is also the desired final state) during the 
second subinterval, by adding the Integral of the square 
of the deviation to the cost during the subinterval. By 
penalizing the deviation from the desired trajectory in this 
way it is hoped that the optimization procedure can be forced 
to choose a switching time, t^, which will give the system 
a rapid rise time and little overshoot. This Cost function 
is plotted as a function of the time to the first switching, 
t^, in Figure l8. Zero Initial conditions were Imposed 
upon x^(t) and X 2 (t) in determining this Cost function'. 

In order to test the usefulness of the function minimi- 
zation approach, the Cost function, J^, was minimized as a 
function of t^ using a i^'oruran IV subx-ouuiiie which implements 
the Tangent Search method for the minimization of a function 
with constraints. See [91 and [10] for further details of 
the program and the method. 

The result of the optimization is shown as Figure 19, 
which shows the "optimum" time response of this second 
order, swltched-parameter system as the functions x^ and 
These response curves should be compared to those of the 
constant-parameter second order system that are shown as 
x^ and x^ which are taken from the response curves of 
Figure The Improvement is evident, as is the price that 

must be paid for "snappy" response; namely, that the time 
response of the velocity (x 2 (t)) becomes impulse-like, indi- 
cating that the system of x^ and x^ is undergoing some 
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Figure l8. Cost Function of the Second Order, Switched' 

Parameter System of Figure 17 as a Function of 
• the Time to the First Switching for a 10 Volt 
Step Input. 
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Figure 19. Time Response of the Second Order, Switched- 

Parameter System of Figure 17 After the Switch- 
ing Time, t]_, has been Optimized. The Optimum 
Value, t]_ = 1.216 seconds. x^_ and X 2 are the 
Responses of the States of the Optimized System. 

and X 2 are the Responses of a Related Constant 
Parameter System. ■ The Input is a 10 Volt Step 
in Both Cases. 
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relatively violent accelerations compared to the soft res- 
ponse of the constant-parameter, linear system. 

The final example demonstrates the Improvement that can 
be made In the step response of the third order, swltched- 
parameter system of Figures 5 and 6 (shown below as Figure 
20) . The time responses of this system to a 10 volt step 
input, for two arbitrary switching times are shown In 
Figure 7 which Is reproduced below as Figure 21. 

Two Cost functions were constructed for this system. 

As for the previously mentioned Cost functions, these two 
were plotted assuming that the system's states started at 
zero Initial conditions. One was dependent on the values 
of the states at the endpoint of the control Interval only, 
analgouc (26) and Figure 9 for the system cf Fimurs 11; 
and the second Cost function was analgous to (27) and 
Figures 9 and 10. The Cost function which depends only on 
the deviations of the endpoint values of the states from 
some desired final value Is shown In Figure 22 for a system 
with a 10 volt step Input. It Is seen to be not markedly 
different from Figure 12 for the 10 volt step curve except 
that the minimum value Is displaced slightly to a larger 
value of t^. The minimum Is seen to be for t^ slightly less 
than 1.6 for an essentially zero value of the Cost function. 
Although It Is not evident from Figure 22, the Cost function 
continues Its rapid rise as t^ Is Increased beyond t^ = 2.0 
in the same manner as that shown In Figure 10, implying that 
the minimum at t^ = 1.6 Is a global minimum. 
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te [0 ,t^] 




Figure 20. Replicas of Figures 5 and 6. Block Diagrams 
of a Third Order, Switched-Parameter System, 
(a) Prior to Switching. (b) After Switching. 
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Figure 21. Replica of Figure 7. Time Response of the Third 
Order, Switched-Parameter System of Figure 20 to 
a 10 Volt Input Step. (a) Switching Time, 
t]^ = 1.699 Seconds. (b) t^ = 1.3^9 Seconds. 
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Figure 22, Cost Function for the Third Order, Switched- 

Parameter System of Figure 20, as a Function of 
the Time to the First Switching for a 10 Volt 
Step Input; VJhere only the Square of the Devia- 
tion of the States from the Desired Final Value 
is Counted as Cost. 
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This Cost function was minimized using the same method 
as was used for the second order system Just considered with 
the addition that three different values of step input 
amplitude were used. The "optimum" switching time was found 
to be the same for all of the inputs, t^ = 1.562 seconds. 

The time responses of the output states, X 2 (t), for the 
optimum response to the various step inputs are shown in 
Figure 23* Figure 24 shows a similar plot of the time 
responses of x^(t), the derivatives of the output states. 

It should be noted from Figure 23 that the character of 
the optimum response is independent of the input step 
amplitude, the various responses differing only in magni- 
tude. This result is not unexpected. It is a reflection 



/lie o.yoL/vJi 






independence of t^ of the values of the states during the 
control interval. Just as the third order system considered 
at the beginning of this appendix showed linear behavior 
for the same reasons. 

In order to obtain some feeling for the influence the 
choice of Cost function has on the "optimum" response 
obtained, or, what is the same thing, on the value of t^ 
produced by the minimization of the Cost function, another 
Cost function was constructed for the same system. This 
new function was the analog of (27) and Figures 9 and 10. 
The new Cost function is plotted as a function of the time 
to the first sv/ltching in Figure 25* There, it may be 
noted that the "optimum" value of t^ is not the same as 
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Figure 23* Optimum Time Response, $ 2 (t), of the Third Order, 
Swit ched-Parameter System of Figure 20 Obtained . 
by Minimizing the Cost Function of Figure 22 
with Respect to the Time to the First Switching, 
for Various Amplitude Input Steps. Optimum 
Switching Time, 't^ = 1.562 Seconds. 










f- 








Figure 24. 



Optimum Time Response, x^(t), for the Same 
System as in Figure 23- 



85 







I 




COST 





Figure 25- Cost Function for the Third Order, Switched- 

Parameter System of Figure 20 as a Function of 
the Time to the First Switching for a 10 Volt 
Input Step, where the Square of Both the Devia- 
tion of the States from the Desired Endpoint 
Values and the Deviation of the States from the 
Desired Final State During the Second Subinter- 
val are Counted as Cost. 
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it was for Figure 22, but Is increased in value slightly to 
about 1.6 seconds. 

Again, a function minimization was done on the Cost 
function. The "optimum" switching time so obtained was, 
t^ = 1.600 seconds. This optimization was also accomplished 
for the same three amplitudes of input step as were used 
in Figures 23 and 24. The "optimum" trajectories, x^Ct) 
and x^Ct), are plotted in Figures 26 and 27. 

The main difference to be noted between the responses 
of Figures 23 and 26 is that the first "optimum" traj’ectory 
shows a slight undershoot, and the second a slight over- 
shoot . Both optimum systems show a marked Improvement over 
that shown by the same system in Figure 21. As was the case 
fnr r.he 1 i rr.r. Cost furi'ccion consiat;red i oi‘ Unxo svstsm, ohc 
optimum switching time is seen to be Independent of the 
input step function amplitude, and for the same reasons. 

Cost function minimization is, thus, seen to provide a 
useful means of "optimizing" this particular kind of 
sv;ltched-parameter system. The utility of this kind of 
procedure for Improving the responses of other kinds of 
switched parameter systems remains to be demonstrated. 
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Figure 26. Optimum Time Response, Xp(t), of the Third 

Order, Swlt ched-Parameter System of Figure 20 
Obtained by Minimizing the Cost Function of 
Figure 25 with Respect tp the Time to the 
First Switching, for Various Amplitude Input 
Steps. Optimum Switching Time, = I. 6 OO 
Seconds . 
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Figure 2?. Optimum Time Response, x-(t)j for the Same 
System as In Figure 26. 
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